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ABSTRACT 

We propose magnetic proto-neutron star (PNS) winds driven by Alfven waves as well as the neutrino 
heating as an appropriate site for the r-process nucleosynthesis. Alfven waves excited by surface motions 
of a PNS propagate outwardly, and they heat and accelerate the wind by dissipation. Compared with the 
wind purely driven by the neutrino heating, larger entropy per baryon and shorter dynamical time scale 
are achieved, which favors the r-process. We study reasonable cases that a wave amplitude is 10% of 
Alfven speed at the surface to find that a PNS with surface field strength, > 5 x lO^^G, gives suitable 
wind properties for the r-process, provided that dissipation length of the wave is at most ~ 10 times of 
the PNS radius. We also compare properties of transcritical and subcritical winds in light of the r-process. 
We finally discuss possibilities of detections of 7-rays from radioactive nuclei and absorption lines due to 
Ba in supernova remnants which possess magnetars. 

Subject headings: magnetic field - nuclear reactions, nucleosynthesis, abundances - plasmas - stars:neutron - super- 
novae: general - waves 



1. Introduction 

An origin of elements produced via rapid neutron 
capture process (r-process) still remains a mystery. It 
is one of the most important astrophysical problems to 
be solved. First, most massive nuclei (up to mass num- 
ber, A ~ 250) in the universe are synthesized through 
the r-process. It is not only scientificly interesting but 
also important to explore how, where, and when they 
are synthesized. Second, some r-process nuclei can be 
used as chronometers (Cowan et al. 1999; Wanajo et 
al. 2003). For example, the half-Uves of ■^■^^Th and 
^'^^U are 1.405 xlO^o yr and 4.468x10^ yr, respec- 
tively. We can estimate the ages of metal-poor as- 
trophysical objects by their observed abundance ratio 
provided that the mass-spectrum of these produced r- 
process nuclei is predicted precisely. Third, chemical 
evolution of the r-process elements gives us a clue to 
reveal the history of the evolution of our Galaxy itself 
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(Ishimaru & Wanajo 1999; Tsujimoto, Shigeyama, & 
Yoshii 1999; Argastetal. 2004). 

The conditions for the successful r-process nucle- 
osynthesis are (Hoffman, Woosley, & Qian 1997; here- 
after HWQ97): (i) high density, (ii) high entropy per 
baryon, (iii) short dynamical timescale, and (iv) small 
electron fraction, Y^. This is because r-process nuclei 
are synthesized through the non-equilibrium process 
of the rapid neutron capture on seed nuclei, products 
as a consequence of the a-rich freezeout (HWQ97). In 
other words, an explosive, and neutron-rich site with 
high entropy is a promising candidate for the r-process 
nucleosynthesis. 

The most probable sites, which possibly meet 
the above conditions, are collapse-driven supernovae 
(SNe; e.g., Woosley et al. 1994) and/or neutron star 
(NS) mergers (Freiburghaus et al. 1999; Rosswog 
et al. 2000). However, the collapse-driven SNe 
are thought to be more probable sites, because metal 
poor stars have akeady contained the r-process nuclei. 
Mcwilliam et al. (1995) reported that Eu, one of the 
r-process elements, is detected in 1 1 out of 33 metal- 
poor stars, which proves that r-process nuclei have 
been already produced at the very early stage in our 
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Galaxy. Comparing the event rate of collapse-driven 
SNe (10~^ yr~^ Gal~^; van den Bergh & Tammann 
1991) with that of the NS merger (IQ-^ yr"! Gar^ 
van den Heuvel & Lorimer 1996; Bethe & Brown 
1998), we can see that collapse-driven SNe are fa- 
vored since they can supply the r-process nuclei at 
early epochs. Moreover, Cowan et al (1999) reported 
that the abundance ratio of r-process nuclei in metal 
poor stars are very similar to that in the solar system. 
This indicates that the r-process nuclei are synthesized 
through the similar conditions; most of the r-process 
nuclei are from a single candidate. So we assume in 
this paper that most of the r-process nuclei are synthe- 
sized in the collapse-driven SNe. 

Many works have been done to uncover the nature 
of the r-process nucleosynthesis in the collapse-driven 
SNe so far. However, it seems quite difficult to pro- 
duce the r-process elements compatible with the ob- 
servations. Calculations of the neutrino-driven winds 
under the Newtonian gravity show that entropy per 
baryon is too small to succeed in the r-process (Taka- 
hashi et al.l994; Qian & Woosley 1996; hereafter 
QW96). In order to solve this difficulty, QW96 fur- 
ther included a post-Newtonian correction to the grav- 
itational force. Cardall and Fuller (1997) extended 
QW96 by considering a fully general relativistic treat- 
ment. Both showed that the relativistic effects fa- 
vored the r-process nucleosynthesis. After an exten- 
sive survey of the general relativistic effects, Otsuki 
et al. (2000) gave a quantitative conclusion that the 
r-process was realized in the strong neutrino-driven 
winds (Ly ~ lO^^erg s~^) as long as a massive 
(~2.0 Mq) and compact (~ 10 km) proto-NS (PNS) 
is formed. This is finally confirmed by Sumiyoshi et 
al. (2000), Wanajo et al. (2001), and Thompson, Bur- 
rows, & Meyer (2001). Although it is very intriguing 
as the solution is firstly discovered by the introduction 
of the relativistic effects, the above conditions can be 
satisfied only by a soft equation of state (EOS) based 
on a few nonstandard models for sufficiently cold nu- 
clear matter (Wiringa, Fiks, & Fabrocini 1988). In 
fact, the r-process nuclei can not be produced in the 
numerical simulations with a normal EOS (Sumiyoshi 
et al. 2000). Thus it seems that the difficulty can not 
be solved completely by the effects of the general rel- 
ativity only. 

There is only one report of the successful nucle- 
osynthesis (Woosley et al. 1994). In their numeri- 
cal simulation, the entropy per baryon becomes high 
enough to produce the r-process nuclei successfully at 



very late phase of the neutrino-driven wind (^ 10 s 
after the core-collapse). However, such high entropy 
is inconsistent with a reliable analytical formulation 
based on the Newtonian gravity (QW96). The gen- 
eral relativistic effects, taken into account in Woosley 
et al. (1994), cannot explain the discrepancy, ei- 
ther (Otsuki et al. 2000). Furthermore, nuclei with 
A ~ 90 are overproduced, which is incompatible with 
the solar system abundances. The successful mass- 
spectrum at the late phase would also be destroyed by 
neutral-current neutrino spallations of nucleons from 
*He (Meyer 1995). Although the result of Woosley et 
al. (1994) is remarkable, it would not be concluded 
that the problem is settled completely. 

Facing such difficulties, it is worth exploring addi- 
tional effects that encourage the r-process nucleosyn- 
thesis. It is argued that the wind properties are im- 
proved to be suitable for the r-process by a small in- 
crease of energy input at a few tens kilometers above 
the surface, which may be operated by wave mecha- 
nism and/or magnetic activities (QW96; Thompson et 
al.2001). Recent observations show that some frac- 
tion of NSs has relatively strong magnetic field, B > 
lO^^G (e.g. Guseinov, Yazgan, & Ankay 2003). It 
is expected to effect dynamics and energetics of the 
neutrino-driven winds because the magnetic pressure 
becomes comparable to or even larger than the thermal 
pressure (Thompson 2003). As a result, abundances of 
the r-process products would be different in the winds 
from such magnetic PNSs. Let us introduce previous 
works dealing with the effects of magnetic field around 
the PNSs. Nagataki (2001) pointed out importance of 
the magnetic field and investigated a steady, subsonic, 
and rigidly rotating jet, although the r-process condi- 
tions can not be satisfied in the winds. Nagataki and 
Kohri (2001) further examined a steady and subsonic 
wind by using magneto-rotational wind solutions (We- 
ber & Davis 1967), although they only explored field 
strengths, 5 x 10^^ G, which is not dynamically impor- 
tant because of the complicated critical point topolo- 
gies encountered in the steady state solution. As a re- 
sult, they found that the wind properties are still inap- 
propriate for the r-process nucleosynthesis. Thompson 
(2003) considered an ordered dipole with larger field 
strength up to lO^^G and pointed out an importance of 
trapping of the plasma heated by the neutrino. He esti- 
mated enhancement of the entropy due to the trapping 
to conclude that the ordered dipoles with surface fields 
of > 6 X 10^"'G give the suitable circumstances for the 
r-process nucleosynthesis. 
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We propose another mechanism operated by mag- 
netic fields. Alfven waves are supposed to play a 
role in heating and accelerating astrophysical plasma. 
In particular those generated by surface convective 
motions are reliable sources in the solar corona and 
wind. As the waves have momentum and energy, 
their dissipation leads to acceleration and heating of 
the surrounding plasma (e.g. Belcher 1971; Suzuki 
2004). Excitation of Alfven waves is also introduced 
in NSs with strong magnetic field (Thompson & Dun- 
can 1996). Then, those Alfven waves will contribute 
to the heating and acceleration in the magnetosphere, 
similarly to the solar corona. To date, however, a role 
of these waves in PNS winds has not been studied. 
Then, we investigate effects of the Alfven waves on 
the physical conditions such as the entropy per baryon 
and the dynamical timescale, and we discuss whether 
they can help the synthesis of r-process nuclei. 

In section 2, we explain the formulation for Alfven wave- 
driven PNS winds in the hot bubble. In section 3, we 
show the results. Discussions and Conclusion are pre- 
sented in section 4 and 5. 

2. Methods 

2.1. Formulations 

We present formula for the PNS winds in existence 
of the Alfven waves which excited at the surface. Al- 
though this is not addressed in this paper, excitation 
of Alfven waves is studied in NSs with strong mag- 
netic field by several authors (Thompson & Duncan 
1996; Thompson, Lyutikov, & Kulkarni 2002; Ko- 
jima & Okita 2004). A certain fraction of magnetic 
energy in the core of a NS is transferred to seismic 
waves in the crust. Such waves will be coupled with 
Alfven waves in the magnetosphere (Blaes et al. 1989). 
In other words, Alfven waves would be excited from 
around the surface. Moreover, the PNS is convectively 
unstable(Epstein 1979). Drived turbulent- like motions 
(Keil, Janka, & Muller 1996) will further contribute to 
the excitation of the Alfven waves. 

We consider one-dimensional magnetic flux tube 
which is radially open. Alfven waves excited with 
amplitude of magnetic filed, SBq, at the PNS surface 
propagate outwardly along with this field fine. The 
plasma also flows along it when the ideal MHD con- 
dition is satisfied. Then, it can be assumed that all the 
physical variables depend only on radial distance, r. 
We also assume steady state with continuous injection 
of the Alfven waves from the surface. We do not con- 



sider effects of the rotation of a PNS and the general 
relativity to focus on a role of Alfven waves in this 
paper. 

Conservation of magnetic flux, V • B = 0, gives 



Br^ = Bgrg = const, 



(1) 



where B is radial field strength and Bq is field strength 
at the PNS surface, r = ro. Mass conservation under 
the steady state condition becomes 



const, 



(2) 



where p is density and v is velocity. Momentum con- 
servation is expressed as 



dv 
dr 



GM 



1 dP 
p dr 



1 dR,: 



p dr 



(3) 



where M is mass of the PNS and P is total pressure of 

nonrelativistic matter, relativistic particles, and photon 
radiation (QW96). P^ is wave pressure which can be 
written as a function of wave amplitude in magnetic 
filed, SB, as P^, = ^ (Lamers & CassineUi 1999). 
An equation for internal energy, £,is 



,d£ P dp. 

Qv + Qw =v{- ^—), 

p^ dr 



~ dr 



(4) 



where denotes neutrino heating and gi„ is heating 
due to dissipation of the Alfven waves. 

In this study, we consider three processes with re- 
spect to neutrino heating and cooling, that is, neutrino 
scattering processes on electrons and positrons, neu- 
trino absorption on free nucleons, and electron and 
positron capture on free nucleons. can be calculated 
as a function of neutrino luminosity, Ly, and energy, 
tv, following Tubbs & Schramm (1975), Goodman, 
Dar, & Nussinov (1987), and QW96. We do not con- 
sider heating from neutrino/anti-neutrino annihilation 
and cooling from electron/position annihilation which 
are taken into account in QW96, because they are not 
effective as shown by Otsuki et al. (2000). We set (7,^=0 
for T <0.5 MeV, because free nucleons are bound into 
a-particles and heavier nuclei and electron-positron 
pairs annihilate into photons. 

The wave heating, qw{> 0), satisfies an equation 
describing variation of wave energy: 



-P(lu 



dP^ 
dr 



va 



Va + V 



(5) 



where va = B/^/A-Kp is Alfven speed and = 

Stt 



^^{va + ^v) is wave energy flux (Jacques 1977). 
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is wave action constant which is expresses as 
(Jacques 1977) 



_ 5B^ (VA + V){VA + V) 



(6) 



which is a conserved quantity in moving media if the 
wave does not dissipate. Physical meaning of eq.(5) is 
quite similar to the 1st law of thermodynamics; vari- 
ation of wave energy flux, V • F^, is determined by 
work done by wave, —v^^, and wave dissipation, 
—pQw Alfven wave pressure gradient is written as a 
function of q'w as (Lamers & Cassinelli 1999) 



1 dR,: 



5B^ 3w + WA 1 dp 



p dr 



2{v + Va) 327r/9 v + va p dr 



(7) 



Alfven waves are supposed to hardly dissipate due 
to the uncompressional character, if they travel in uni- 
form media. However, amplitude of those propagating 
in the stratified plasma, such as atmosphere of stel- 
lar objects, is amplified so that non-linear damping 
processes becomes important. Ingoing Alfven waves 
are excited by parametric decay instability, coupling 
between the outgoing waves and density fluctuations 
(Goldstein 1978). These ingoing Alfven waves in- 
teract with the pre-existing outgoing waves to lead 
to turbulent-like cascade and subsequent wave dissi- 
pation. The parametric decay also generates sound 
waves, which easily dissipate due to their compressive 
character. Linearly polarized waves dissipate through 
formation of fast MHD shocks in magnetically dom- 
inated plasma as a consequence of nonlinear steepen- 
ing of flie wave fronts (HoUweg 1982; Suzuki 2004). 
Transverse gradient of Alfven speed leads to resis- 
tive and viscous dissipation due to phase mixing of 
Alfven waves traveling along neighboring field lines 
(Heyvaerts & Priest 1983). Kinetic effects, such as 
ioncyclotron resonance, might work in the dissipation 
(e.g. Craimier, Field, & Kohl 1999). Observation in 
the solar corona shows an evidence of dissipation of 
Alfven waves in a region closed to the surface (Doyle, 
Teriaca, & Banerjee 1999). Therefore, Alfven waves, 
if they exist, are also expected to dissipate in the PNS 
winds, though the dissipation mechanisms are quite 
complicated as discussed above. In this paper, we 
adopt a simple parameterization for the wave damping 
by using a dissipation length, I, 



,oexp 



(8) 



where H^fl is wave action at the PNS surface, r ~ tq. 
Equations (5) and (8) give an exphcit form of heating 
function as 

.IVA rl exp(^) 

(Iw — ; J^w,0-n ; (y) 

pv + VA ?■ t 

Once initial wave energy flux, F^ o(— ^fui.o for w <C 
Va), is fixed, / controls distribution of energy and mo- 
mentum transferred from the waves; heating and accel- 
eration occurs in inner regions for cases with smaller 
I, and vice versa. 

Let us give a rough estimate of I based on analogy 
to the nonlinear damping of Alfven waves in the solar 
corona. Generally, / is expected to have a linear de- 
pendence on wave length, A = vat, where r is wave 
period: 

l = fX, 

where / is a proportional coefficient which is ex- 
pected to have a negative dependence on 5B/B for 
the nonlinear damping. Suzuki (2004) considered 
Alfven wave with A ~ lO^km (va ~ lOOOkm/s and 
T ^ lOOs) as a typical one in the solar corona. He 
studied dissipation by fast shocks, one of the nonUn- 
ear mechanisms, to find that I becomes an order of Rq 
for SB/B ~ 0.2 - 0.3, where Rq = 7 x lO^km is 
the solar radius. This shows / 10, which means 
that the Alfven waves dissipate typically by 10 times 
of A. In PNS winds, timescale for convections can be 
regarded as one of the typical wave periods. Hydro- 
dynamical simulations show that convective cells with 
several kilometers move with velocity 10^ — 10''km/s 
(Keil et al.l996), which gives r ~ 1ms. We consider 
Bo > lO^^G for Alfven wave-driven winds, which 
corresponds to va lO^km/s. Then, typical wave 
length becomes A = vat ~ 10km(= ro). If the am- 
plitude is similar to that in the solar wind, the sim- 
ilar /(~ 10) may arise. It follows that we can as- 
sume I ~ lOro as a fiducial value. However, / would 
vary a lot if waves with different r dominantly worked. 
Therefore, we consider three cases, I = (5, 10, 30) x ro 
in this paper. 

2.2. Numerical Integration 

We describe a practical method to determine the 
wind structures for various inputs of Alfven waves. 
Adopted parameters for a PNS and the input neutri- 
nos are summarized in tab.l. We set the density at 
the iimer boundary to be 10^° g/cm^, following the re- 
sult of Wilson's numerical simulation in Woosley et 
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al. (1994) (see also Otsuki et al. 2000). Temperature 
at the surface is derived from energy balance between 
the neutrino heating and cooling, following QW96: 



To = 1.19 X 10 



10 



D 3,3 

,51 -"-6 ^Pe.MeV 



[K], 



(10) 



where Li^^^i is the individual neutrino luminosity in 
10^^ ergs s~^, Rq is the neutron star radius in 10^ cm 
(= 10km), A = 1.293 MeV is the neutron-proton mass 
difference, and e,y mcV is a neutrino energy in MeV. 
We assume that the neutron star radius is equal to the 
neutrino sphere radius. 

Amplitude of Alfven waves is set to be SBq/Bq = 
0.1 at the surface, although it is quite uncertain. 
SBq/Bq = 0.1 indicates that velocity amplitude is 
10% of the phase speed (Alfven speed in this case), 
which is often seen in astrophysical objects; for ex- 
ample, at the solar surface speed of the convective 
motions (~ Ikm/s) is ~ 10% of the sound speed 
(^ lOkm/s). By the assumption of the constant 
SBq/Bq the input wave energy flux, Fyj^ ~ -ffw,o. 
is proportional to Bq (note that va,o oc Bq). 

Integration starts from the surface for an initial 
guess of velocity, vq. Temperature is integrated ac- 
cording to eq.(4). Background magnetic field, B, is 
fixed by the conservation of magnetic flux (eq.(l)). SB 
is determined in order to satisfy the dissipation model 
of Alfven wave (eq.(8)). We integrate velocity by an 
equation, transformed from the momentum equation 
(3) with help of eq.(7), 

1 dv 
V dr 



^ r ' 3v 



+ 



2(v + vk) 
(11) 



where = (f-^)^^^ is adiabatic sound speed and 



w 



I dB^ 3v + VA 

32'Kp V + VA 



is a variable in unit of velocity, indicating contribu- 
tion of the Alfven waves. When deriving eq.(l 1) from 
eq.(3), we have assumed that pressure by relativis- 
tic electrons and positrons and photon radiation domi- 
nates the other components, and the degeneracy is neg- 
hgible, namely, 

P^^^, (12) 
4 3' ^ 



where a(= 7.56 x 10^^"''erg cm^'^K^'*) is a Stefan 
Boltzmann radiation constant. The above assumption 
is applicable only for T > 0.5MeV because otherwise 
the electron-positron pairs annihilate and the remained 
electrons become nonrelativistic. Although the effect 
of the pair annihilation gives a modification of eq.(12) 
by a factor of (11/4)^/''', this seems to influence on the 
dynamics quite httle (Sumiyoshi et al. (2000)). Con- 
tribution of the non-relativistic gas to the pressure is 
also small for T > O.lMeV in our Alfven wave-driven 
winds. This is because density in the Alfven wave- 
driven winds is smaller than the neutrino driven-winds 
on account of rapid acceleration. Once v is determined 
by eq.(ll), p is derived from the mass conservation 
equation (2). 

For the velocity structure we both consider trans- 
critical and subcritical solutions. When deriving trans- 
critical solutions, a critical point is derived from the 
condition that both numerator and denominator of 
eq.(ll) are zero. Then, we iteratively determine an 
initial guess of vq to smoothly pass through the crit- 
ical point by a shooting method. On the other hand, 
subcritical solutions cannot be uniquely determined. 
We derive subcritical solutions for given mass flux by 
reducing vq from the transcritical values. 

We have checked an accuracy of the integration by 
monitoring whether how precisely a Bernoulli equa- 
tion shown below is fulfilled as a function of r: 



\P TS_ 

2 jTiN r 



GM Fw 

+ — 

pv 



dr — , 



(13) 



where is the nucleon rest mass and S = ton ^ is 
entropy per baryon (for non-degenerate and relativistic 
gas). We have found that our numerical integration 
satisfies the above equation at least within 3% which is 
acceptable for the purpose of our model calculations. 



Table 1 : Adopted Parameters, ro is PNS radius, po is 
surface density, and M is PNS mass in unit of solar 
mass, Mq. Ly is neutrino luminosity per flavor and 
is neutrino energy. 





10km 


pa 


10i"gcm-3 


M 


1.4M0 




lO^^erg s-i 




12MeV 




22MeV 




34MeV 



5 



2.3. Condition for Entropy 

When we discuss whether r-process nucleosynthe- 
sis can occur or not, we use a criterion given by 
HWQ97, which is written as 

// \l/3 

5 > 2 X lO^Ye ( ] (14) 

for >0.38, where S is measured in unit of Boltz- 
mann constant, fcs, and texp is expansion time scale 
defined as time elapsed from T = 9xlO''KtoT = 
2.5 X lO^K. This is the condition for production of the 
r-process nuclei with mass number A -^200 (QW96). 

Electron fraction Yg at the time when a-rich freeze- 
out takes place can be estimated as (QW96) 

y _ A _|_ Lp^ ePe,McV — 2A + 1.2A^/ep_j3icv\ 
''"V e,.e,McV + 2A -t- 1.2A2/e^,,McV / 

When we substitute the parameters shown in tab.l, 
we estimate Yg = 0.43, at which the r-process cri- 
terion (eq.(14)) is appUcable. In reaUstic situations, 
L,y and e,y may vary, hence, we consider cases with 
Ye = 0.4 - 0.5 

3. Results 

In this section, we examine how the dissipation of 
the Alfven waves modifies the wind structures. In par- 
ticular we discuss possibilities of the r-process nucle- 
osynthesis in light of the criterion, eq.(14), for S and 

^exp' 

The transcritical solutions are generally stable if 
sufficient energy injection is given and the winds blow 
into free space. In the PNS winds, however, ejected 
shells exist in the outer region and interaction between 
the shells and winds leads to formation of termina- 
tion shock. We expect the flow stays transcritical if 
a position of the shock is far outside of the critical 
point. Otherwise, the shock works as an outer bound- 
ary so that the flow might become subcritical eventu- 
ally. Thus, we firstly study properties of transcritical 
flows in detail (§3.1), and then, show different aspects 
of subcritical flows for comparison (§3.2). 

3.1. Transcritical Flows 

3.1.1. Effects of Alfven Waves 

Figure 1 illustrates how the wind properties change 
by the effects of Alfven waves. As for Alfven wave. 



we consider two cases of _Bo = 3 x 10 ^^G (dotted) and 
Bq = 5x 10 ^^G (solid). In both case the same dissipa- 
tion length, I = lOro, is adopted. Wind purely driven 
by the neutrino heating is also shown for comparison 
(dashed). Since F^^o oc Bq, the input Alfven wave 
energy flux of case with Bq = 5 x 10^^ is larger by 
a factor of 5 than that of case with i?o = 3 x 10^'*. 
Wind structures (p, T, & v), entropy per baryon (S), 
and heating (g^, & Qyj) are compared in top, middle, 
and bottom panels, respectively. 

One sees in bottom panel that the wave heating is 
distributed around ~ l{= lOro), in contrast to the neu- 
trino heating which is localized close to the surface. 
Then, the wind structures shown in top panel are mod- 
ified in relatively outer region by the Alfven waves. 
A clear difference is seen in velocity structure. The 
Alfven waves directly accelerate the materials by wave 
pressure gradient and the faster winds are attained as 
input energy of Alfven wave increases. Accordingly, 
iexp = 8.9 X 10~^s in the case with larger input of 
Alfven waves is much shorter than texp — 4.4 x 10~'^s 
in the case without the waves. This indicates that 
the wave dissipation reaUzes non-equilibrium circum- 
stances for the nuclear reactions, which encourages 
the r-process nucleosynthesis. Density shows more 
rapid decrease by the presence of the Alfven waves 
to satisfy the mass conservation (eq.(2)). Temperature 
is also slightly lower because adiabatic coohng effec- 
tively works due to the rapid expansion (i.e. accelera- 
tion) of the plasma. 

Middle panel shows that the Alfven waves give 
a great influence on 5* in the wind. First, S is no 
more a constant even in the outer region where the 
neutrino heating is inefficient because the wave heat- 
ing still occurs there. Second, S is larger except in 
a region closed to the surface for larger injection of 
Alfven waves. This trend needs to be interpreted in 
detail because it is opposed to that anticipated from 
the neutrino heating (QW96) . Variation of entropy is 
given by 

dS m^q„ mNfc ,,,, 
-d^ = ^ + ^' ^^^^ 
which is derived from the energy equation (4). For the 
neutrino heating, larger leads to larger v and T in 
the winds. As a result, ||t itself becomes smaller for 
larger and final S is also smaller (QW96). On the 
other hand in the wave heating, an increase of w as a 
consequence of an increase of is smaller than that 
of qw itself as described below. 

Based on our wave heating model, eq.(5), q^ is pro- 
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portional to Bq for the assumed constant initial ampli- 
tude 6Bo/Bq = 0.1. As for v, a rough dimensional 
estimation can be done in a following way. Without 
Alfven waves, velocity, vi, at an arbitrary position 
r = ri can be derived from the BernouilU equation 
(13) as 



2 mN 



TiSi GM 



ri 



GM 
+ 



Jrn V 



(17) 



where subscript '1' denotes physical variables at r\ 
and we have used the fact that the gravitational po- 
tential term dominates at the surface, r = ro- Let us 
assume that velocity increase to vi + Avi by an input 
of Alfven waves with energy flux, F^,o. at ro. Then, 

K + Av,)^ ^ (Ti + Ari)(gi + Agi) GM ^ 

2 rriN ri pivi 



GM 



ar— 
v' 



Fw,o 



(18) 



ro Jro ^ Povo 
where ATi and A 5*1 are modifications of Ti and 5i, 
and we use v' instead of v in the integral of the neutrino 
term because the velocity changes from that in the pure 
neutrino-driven wind (eq.(17)). 

The wave energy converts to both kinetic energy 
(first term on the left of eq.l8) and enthalpy (second 
term) of the flow. We assume that a fraction, e, is 
transferred to the kinetic part, whereas e is automati- 
cally calculated in the numerical integration to satisfy 
the wave energy equation (5). Comparing eqs.(17) and 
(18), we have 



2 2 



w.l 



PoVq piVi 



(19) 



where we neglect the difference of the velocity appear- 
ing in the neutrino term. If Avi <C t^i, ^ in eq.(16) 
simply proportional to oc Bq, hence, an increase 
of Bq raises S. If Avi ^ vi, equation (19) gives 
At;i oc oc By^ for constant vo, provided that 

^VO ^ jjjg j.g^J sj];Qjj];iQ[j y Jj^g ^ pOSltlVC 

povo plvi " r 

correlation with Bo so that the dependence of Avi on 
Bq is weaker. Therefore, ^ in eq.(16) is at least in 

proportion to B^^ , showing that larger field strength 
directly leads to an increase of S. In the Alfven wave- 
driven winds, temperature is also sUghtly lower due 
to the larger adiabatic cooling (fig.l). This effect fur- 
ther contributes to the increase of S as expected from 
eq.(16). In conclusion, Alfven waves effectively raise 



entropy in the winds to give the suitable circumstances 
for the r-process. 

Since S varies in a region where the r-process oc- 
curs, it is not straightforwards whether the r-process 
criterion, eq.(14), is directly applicable. In this pa- 
per, we adopt 5 at T = 0.2MeV as a necessary con- 
dition, based on the fact that a ratio of neutrons to 
seeds, a key to determine the final abundance of the 
r-process elements, is once fixed through the a-rich 
freeze-out in T > 0.2MeV. The figure exhibits that 
S{T = 0.2MeV) (filled triangles) increases for larger 
input of the Alfven waves. In particular, case with 
Bo = 5 X lO^-^G gives S = 170, which satisfies the 
criterion, eq.(14), for the r-process nucleosynthesis, if 
Ye < 0.41. 

3.1.2. Dependence on Dissipation Length 

Figure 2 compares the wind properties for differ- 
ent dissipation lengths, Z = 5 & 30 but the same field 
strength, Bq = 5 x lO^^G. As illustrated in bottom 
panel, I controls the distribution of heating; the waves 
dissipate more efficiently around / so that the heat- 
ing and acceleration (not shown) are most effective 
there. For smaller I, the Alfven waves influence on the 
wind in an inner part; higher v, smaller p, and T are 
obtained as seen in top panel, which is consistent with 
the trends obtained when considering larger wave heat- 
ing ( § 3 . 1 . 1 ). 5 is also larger in the inner region (middle 
panel) as the wave heating is enhanced there. How- 
ever, these trends become reverse in the outer region 
as the wave dissipation becomes efficient for larger I 
case. 

The wind properties in T > 0.2MeV are important 
with respect to the r-process nucleosynthesis, since the 
neutron-to-seeds ratio is determined there. Wind struc- 
ture is greatly improved there in case with Z = 5 to 
give smaller ioxp and larger S for the successful nucle- 
osynthesis, while the modification is smaller in case 
with Z = 30 as the wave dissipation occurs in farther 
outside. Therefore, rapid wave dissipation (short I) is 
favored for the r-process. 

3.1.3. Bq & I for R-process 

In fig. 3 we show result of the wind properties for 
various Bo and I in texp-S{T = 0.2McV) plane, over- 
layed with the r-process criterion (eq.(14)) for = 
0.4 and 0.5. We do not show results with Bq > 6 x 
lO^^G, because both wind velocity and Alfven speed 
are approaching to the Ught speed, c, and our non- 



7 



relativistic treatment is not valid. (The final flow speed 
becomes ~ (0.4 — 0.5)c and the Alfven speed at 
r ~ 5ro becomes ~ c for Bq = 6 x lO^^G.) The 
figure shows that the Alfven waves play a role for 
-^0 ^ 10^'*G. Larger Bq and smaller I lead to larger 
S{T = 0.2MeV) and smaller texp as described in 
§3.1.1 and §3.1.2. Models with Bq > 5 x 10"G and 
I < lOro satisfies the r-process condition, eq.(14), for 
Ye = 0.4. The condition for = 0.5 is also satisfied 
with sHghtly larger Bq {= 6 x lO^^G) if I < lOro, as 
dependence on Bq is sensitive. 

The value, Bq — 5 x lO^^G, can be understood by 
a simple energetics consideration. It is useful to intro- 
duce a plasma /3 value, (3 = P/{B'^ /Sw), a ratio of 
plasma pressure to magnetic pressure. Let us consider 
/3atr = lOro: 

. ,1 f Bo V /0.15MeV\'' 

where subscript '10' denotes quantities measured at 
r = lOro. Here, we have used the conservation of 
the magnetic flux (eq.(l)) to relate Bio with Bo, and 
we assume that the radiation pressure dominates in P. 
Equation (20) shows that the flow at lOro is domi- 
nated by the magnetic field energy by 10 times for 
Bo = 5x lO^^G. Transfer of small fraction of the mag- 
netic energy to the radiation plasma gives a great in- 
fluence on the energetics and dynamics. Alfven waves 
play a role in this energy transfer in our framework, 
and the wind properties are improved to be suitable for 
the r-process nucleosynthesis. Sensitive behaviors of 
the wind parameters on Bo seen in fig. 3 can also be 
interpreted by squared dependence on Bo in eq.(20). 

3.2. Subcritical Flows 

We have discussed the r-process nucleosynthesis in 
the transcritical flows so far This is reasonable if the 
plasma expands into space with much lower density. 
In the real situations, however, the PNS winds would 
be interrupted by the ejected shells existing in the outer 
region to form the termination shocks. If a position of 
the shocks are not far from the critical point, the winds 
might become subcritical flows eventually. Terasawa 
et al. (2002) pointed out that these dynamics at the 
outer boundary gave an influence on the r-process nu- 
cleosynthesis. Hence, it is worth studying the dynam- 
ics of the subcritical flows with Alfven waves. In this 
subsection, we examine how the structures in the sub- 
critical flows are modified compared with the transcrit- 



ical case and argue possible influences on the r-process 
nucleosynthesis. 

Infinite numbers of subcritical solutions exist for 
different choices of vo- Physically, a certain solution 
is selected by conditions at the outer region, such as 
pressure in the shell ejected by the SN explosion. The 
shell moves outward and the pressure would decrease 
as it expands. Then, it is required to treat the outer 
boundary in a time-dependent manner (e.g. Sumiyoshi 
et al. 2000) for accurate calculations. However, we 
construct subcritical wind structures under the steady 
state condition by reducing vo in an ad hoc way for 
simplicity's sake. 

Figure 4 compares subcritical flows with the tran- 
scritical one for Bo = 5 x and I = lOro. As 
for the subcritical flows, we show cases with smaller 
Vo by 0.1%(case A), 1%(B), and 2.5%(C) than that of 
the transcritical case. The figure shows that the wind 
structures in the iimer region (r < lOro) are almost 
identical, so that t^xp differs only by a factor of 1.6 
even between case C and the transcritical flow. Clear 
difference is seen in the outer region. Speed of subcrit- 
ical flows in the outer region is slower, by definition, 
than the transcritical flow. Accordingly, the density is 
much higher to satisfy the mass conservation. Tem- 
perature is also higher because the wave dissipation 
works in the plasma heating rather than the wind ac- 
celeration. The difference in the entropy (oc T^/p) is 
small because the change of p is compensated by that 
of T, whereas the subcritical cases give slightly larger 
S{T = 0.2MeV). The r-process criterion, eq.(14), is 
influenced quite little as the tiny variances of S and 
<exp further cancel out. 

We would like to discuss implications on the r- 
process nucleosynthesis inferred from these different 
wind properties. The r-process nuclei are synthesized 
by following two sequential processes (Woosley & 
Hoffman 1992). In the inner region with higher tem- 
perature (T > 0.2MeV), a-process proceeds to syn- 
thesize the seed nuclei with mass number, A ~ 100. 
In the outer region with lower temperature (T < 
0.2MeV), reactions concerning charged particles cease 
and the neutron capture by the seed nuclei, namely 
the r-process in the narrow sense, dominantly occurs 
to produce the heavy nuclei. Since the dynamics in 
the inner region resemble each other, nucleosynthesis 
by the a-process is expeted to be similar so that the 
neutron-to-seed ratio fixed at T ~ 0.2MeV is not dif- 
ferent so much whether the flows are transcritical or 
subcritical. 
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A main contrast would arise in the final neutron 
capture in the outer region. An advantage of the sub- 
critical flows is that the seed nuclei have the suffi- 
cient time to capture all the ambient neutrons due to 
the slower flow velocity (Yamada 2004). Temperature 
(T ~ O.lMeV) of the subcritical winds is also rea- 
sonable for the neutron capture (Wanajo et al. 2002). 
We expect that the subcritical winds might be more fa- 
vorable for the synthesis of the r-process nuclei with 
A ^ 200, whereas detailed calculation of nuclear re- 
action is desired to give the final conclusion. 

4. Discussions 

4.1. Configurations of Magnetic Fields 

In this paper, we have only considered flux tubes 
radially open as a simplest case. If some portions of 
the surface are covered with closed magnetic fields, 
the open flux tubes expand super-radiaUy (e.g. Gol- 
dreich & Julian 1969 for the dipole case). We can raise 
mainly two effects with respect to such super-radial ex- 
pansions. First effect is that it reduces Alfven speed, 
leading to shorter dissipation length of Alfven waves. 
It may result in larger entropy due to an increase of 
the heating in the inner regions. Second effect is that 
it enhances adiabatic loss to reduce the wind speed 
(Suzuki 2004) and the temperature. Then, texp would 
be changed, whereas detailed investigation is required 
to tell how it is modified. 

If the plasma in the closed region affects the PNS 
winds by magnetic reconnections, situations are quite 
different. The plasma in the closed regions are once 
heated by the wave dissipation as well as the neutrino 
heating and their sudden break-up leads to inpulsive 
heating. This may also give suitable sites for the r- 
process nucleosynthesis as pointed out by Thompson 
(2003). 

4.2. PNS Rotation 

Massive stars are known to rotate rapidly (Tassoul 
1978). Recent numerical calculations also show that 
the iron core is rotating prior to the collapse (Heger, 
Langer, & Woosley 2000; Heger, Woosley, & Spruit 
2004). If the rotation is taken into account, jet-Uke 
explosion possibly takes place (Kotake, Yamada, & 
Sato 2003; Takiwaki et al. 2004). It is shown that 
such a jet-like explosion makes the products of explo- 
sive nucleosynthesis change so much (Nagataki et al. 
1997, 1998; Nagataki 2000; Nagataki et al. 2003). 



Although the r-process is not favored in steady, sub- 
critical, and rigidly rotating jets in which a balance 
is maintained between the centrifugal and magnetic 
forces (Nagataki 2001), there will be still a possibility 
of the successful r-process nucleosynthesis in a jet-Uke 
explosion. In particular, we should investigate dynam- 
ical jets without assuming a steady state. Numerical 
simulations with rotation and magnetic field are highly 
desired to investigate such a possibility. 

There will also be a possibility of the r-process nu- 
cleosynthesis in the winds on the equatorial plane. Na- 
gataki & Kohri (2001) have investigated only steady 
and subcritical wind solutions which use the model of 
Weber & Davis (1967), although they are not relevant 
for the r-process. Otherwise if we investigated dynam- 
ical winds without assuming the steady-state, we may 
find a situation where successful r-process nucleosyn- 
thesis occurs, as pointed in Nagataki & Kohri (2001). 

4.3. Possibilities of Gamma-Ray Observation 

Qian, Vogel, & Wasserburg (1998; hereafter QVW98) 
discussed possibilities of detection of 7-ray emitted 
by the radioactive decay of the synthesized r-process 
nuclei in SN remnants (SNRs). They estimated 
7-ray emissions of several isotopes with half -life, 
Thi ^ lOyr, assuming the solar r-process abundances 
are explained by the accumulation of the Galactic SNe. 
The derived values are an order of ~ 10~'^7cm~^s~^ 
if a SN explodes at a distance of lOkpc, which is be- 
low the sensitivity of International Gamma-Ray Astro- 
physics Laboratory (INTEGRAL), but at tiie level of 
future proposed telescope, ATHENA. 

Our results indicate that magnetic PNSs could pro- 
duce more r-process elements than the usual PNSs 
with Bq ~ lO^^G. We would hke to give a rough 
estimate on the 7-ray emissions from these magnetic 
PNS in the similar way to that adopted in QVW98. 
Subgroups of the NSs caUed anomalous X-ray pulsars 
(AXPs) and soft gamma repeaters (SGRs) are inferred 
to have strong magnetic field, ~ (Duncan 
& Thompson 1992; Vasisht & Gotflielf 1997;Kouve- 
liotou et al. 1998; Zhang & Harding 2000). Alfliough 
it is still premature to carry out statistical analysis on 
number fraction of these magnetic NSs, a rough es- 
timate can be done. Comparison between the nearby 
AXPs and SNRs gives, < 1/50 of SNe have magnetic 
NSs (Guseinov et al.2003). Arons (2003) adopted 
more conservative limit for birth rate of the magne- 
tars, 10"^ - 10-3yr-\ which indicates that 0.03 - 3% 
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of the SNe remain magnetic PNSs for assumed rate of 
SN, 0.03yr-i. 

Let us consider an extreme case that all the r- 
process elements have been synthesized by the mag- 
netic PNSs. If we adopt 1% for the number fraction 
of the magnetic PNSs as a fiducial value, the required 
r-process yield per PNS for the solar abundance is in- 
creased by a factor of 100 against the case in which 
all the PNSs produce the r-process elements equally. 
Therefore, expected intensities of 7-ray lines from sev- 
eral isotopes become ^ 10~^7cm~^s~^ for a distance 
of lOkpc. Two AXPs are detected within ~ 7kpc 
(Guseinov et al.2003). Their field strength is derived 
as ~ 8 X IQi^G (IE 1841-045) and ~ 7 x IO^^q 
(1E2259+586) from the P—P diagram, whereas, these 
values could be larger if an effect of non-uniform ro- 
tation is taken into account (Lugones & Bobaci 2004). 
Unfortunately, estimated ages, > lOOOyr, of the SNRs 
(Vasisht & Gotthelf 1997; Parmar et al. 1998) are 
much longer than the half Uves (~ lOyr) of the quoted 
nuclei in QVW98, so that it is impossible to detect 
the 7-ray by INTEGRAL. However, our estimate im- 
plies that 7-ray in future magnetic PNSs associated 
with SNe in the Galactic disk are detectable by IN- 
TEGRAL, whereas the possibility of the occurrence is 
low due to small birth rate 10^'' — 10^"^yr^^). 

QVW98 further considered long lived radioac- 
tive nuclei, ^^SsN with rhi = 1.14 x lO^yr. Their 
advantage is that 7-ray can be observed for a rela- 
tively longer period in spite of the weaker emissiv- 
ities. Then, we do not have to wait another SN to 
explode. We can estimate the 7 ray emission to be 
~ 2 X 10-*^((5M/5 X 10-5Mo)(7kpc/d)27cm-2s-i 
for IE 1841-045 (Sanbonmatsu & Helfand 1992), and 
~ 5.5 X lQ-^{5M/b X IG-^Mo) (4kpc/Q!)27cm-2s-i 
for 1E2259H-586 (Coe & Jones 1992), where 5M is the 
average amount of mass of a radioactive r-process nu- 
cleus produced in a magnetar (QVW98). Thus it will 
be difficult to detect 7-ray lines from these long lived 
radioactive nuclei by INTEGRAL. However, the 7 
rays from these AXPs might be marginally detectable 
by ATHENA whose sensitivity is ^ 10~^ 7cm~^s~-'^ 
at ~ 100 — 700 keV in the future, whereas a large 
uncertainty of 5M should be reduced by theoretical 
efforts. 

4.4. Absorption Features of Supernova Spectrum 

There are reports that absorption lines of Ba[II] and 
Ba[lll] are detected in SN1987A as well as in other 



Type 11-P supernova in the early phase of the explosion 
(Mazzali et al. 1992; Mazzah & Chugai 1995), which 
may indicate that r-process nucleosynthesis occured 
in SN1987A. Tsujimoto & Shigeyama (2001) esti- 
mated the amount of synthesized Barium in SN1987A 
is about 6x10~^Mq. Thus if a SN that contains 
a magnetar explodes in nearby galaxies and ejected 
more r-process elements, the spectrum should show 
clearer absorption features due to Ba. In this case, Sr, 
^ 10% of which is produced by the r-process (note 
that Sr is mostly produced by the s-process), would be 
observable as well (Tsujimoto & Shigeyama 2001). 

5. Conclusion 

We have investigated possibilities of the r-process 
nucleosynthesis in the PNS winds driven not only by 
the neutrino heating but also by the Alfven waves. 
Those excited by the surface motions propagate up- 
wardly to heat and accelerate the winds through their 
dissipation. As a result, larger S and smaller i^xp can 
be accomplished. We have considered the waves with 
initial amplitude, SBq/Bq = 0.1 and construct 1 di- 
mensional spherical symmetric wind structures under 
the steady state condition. We have shown that field 
strength, Bq > 5 x 10^''G, satisfies the criterion for 
the r-process by HWQ97, if the Alfven waves dissi- 
pate around a distance, < 10 times of the radius, which 
is reasonable for the waves with period of 1ms. We 
also compared transcritical and subcritical flows with 
respect to the r-process. We have found that the final 
neutron captures are affected by the types of the flows, 
while the neutron-to-seed ratio after the a-process is 
expected to be similar. 

Our results imply that the magnetic PNS is a more 
appropriate site for the r-process nucleosynthesis than 
PNSs with usual field strength ^ lO^^G. Based on this 
consideration, we have given an estimation of 7 ray 
emissions by the radioactive nuclei from AXPs (and 
SGRs), relevant candidates for the magnetic NSs. Our 
estimates show that 7 rays from the short-lived nuclei 
could be detected in future SNe associated with the 
magnetic PNSs by the present telescope, INTEGRAL, 
if one happens in spite of the low possibility, while 
those from the long-lived nuclei would be marginally 
detectable in the present AXPs by the future telescope, 
ATHENA. 
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Fig. 1. — Properties of Alfven wave-driven winds for 
Bo = 5 X IQi'^G (solid) and 3 x 10"G (dotted). Dissi- 
pation length is set to be i = lOro for both cases. Re- 
sult of neutrino driven- wind (without Alfven waves) is 
also displayed for comparison (dashed), top : Com- 
parison of density in lO^g cm~^, temperature in lO^K, 
and velocity in lO'^cm s~^. middle : Comparison of 
entropy per baryon in unit of /jb- Triangles indicate 
positions where T = 0.2MeV. bottom Comparison of 
heating (10^°erg/g s) by the neutrino (g^) and wave 
dissipation (qw). 
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Fig. 2. — Comparison of the wind properties for differ- 
ent dissipation lengths, Z = 5 & 30 x ro but the same 
Bo = 5 X lO^^G. Each panel is the same as in fig. 1 but 
the horizontal and vertical scales are different. 
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Fig. 3.— tcxp(s) (X-axis) and 5 at T = 0.2MeV (Y- 
axis) for various parameters of Alfven waves. Open 
circles, filled triangles, and open squares are results 
with ^ = 5, 10, & 30, respectively. Numbers denote 
magnetic field strength at the surface. For example 
6(14) indicates = 6 x lO^^G. Results with the 
same Bq are connected by dashed lines. Solid lines 
are the conditions for the r-process by HWQ97 (eq.l4) 
forFe = 0.4 & 0.5. 
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Fig. 4. — Wind properties of subcritical flows in com- 
parison with transonic flow (sohd) for Bq = 5x lO^^'G 
and I = lOro. Dashed, dot-dashed, and dotted lines 
are cases with reducing vq by 0.1%, 1%, & 2.5%, re- 
spectively. Upper and lower panels are the same as in 
top and middle panels in fig. 1, but the horizontal and 
vertical scales are different. 
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